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Abstract. 

Starting from the Hamiltonian for a dimer which includes all the electronic and 
elcctron-phonon terms consistent with a non-degenerate orbital, by a sequence of 
displacement and squeezing transformation we obtain an effective polaronic Hamil- 
tonian. The rcnormalized electronic interactions differ from the results of semiclas- 
sical or perturbative treatments. The properties of the variationally determined 
ground state of two particles in the orbital are discussed for variable dimer length 
in the adiabatic limit. 



1 Introduction. 

Thirty years ago, the electron-phonon interaction 7^ = dUj/d(Rj — R;), resulting 
from the modulation of the hopping amplitude Uj due to the displacement Ui of the 
i-th site from its equilibrium position R^ in the lattice, was indicated as responsible 
of the BCS superconductivity in transition metalsjl). Later on its relevance also to 
the physics of quasi- unidimensional conducting polymers Q and layer compounds || 
was recognized. Interest in its effects, however, was rather sporadic until the seminal 
paper on dimcrisation in unidimensional systems [Q. Since then, the Hamiltonian 

H S SH = "/ij( C L C ja + c ]a C ta){ u j ~ "i) W 

is known as the Su-Schrieffer-Heeger (SSH) model[^|. The suggestion was also ad- 
vanced that it had to do with the bipolaronic superconductivity^, a very exotic 
possibility at that time. Its effect in transition metals was reexamined in Ref. 
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All those issues became hot subjects of debate after the discovery of the high- 
temperature superconductors, promoting many investigations Q on its possible ef- 
fects in the cuprates. However, due to the impossibility of an exact solution on an 
infinite lattice, the study of Hssh has generally been limited to the quasi-classical 
approximation, treating (v,j — itj) as a c- number, and resulting in a modification 
of the hopping amplitude tij linearly depending on 7^: t*j = tij — 7^ (uj — u{). 

A second- order perturbative treatment has been recently published^ in the con- 
text of the equivalence between spinless SSH interaction term and the X- Y spin- 
Peierls Hamiltonian, suggesting a quadratic dependence of t*j on 7.^ . On small 
size clusters the problem allows in principle for an explicit analytical solution with 
(v,j — Uj) properly quantized. Here we shall present a non-perturbative study of a 
two-site model system (a dimer) where the assumed set of electron-phonon interac- 
tions generalizes Hssh by including other terms consistent with one electron in a 
non-degenerate orbital, as discussed in Ref.10. 



2 The model. 

Let us consider a general two-site electron-phonon Hamiltonian H = H e i + H. 



H. 



el—ph 



where, in standard notation for a non-degenerate orbitalfl] 



ph 



H e i = e^(n la +n2a)+^[t+X(n 1 ^ a + n2-a)}(c\ a C2a + H.c) + U(n 1 ^n ll +n 2 ^n2i) 



+ (V - J/2)mn 2 - 2J 



SIS Z 2+-(S+S2 +H.c.) 



+ P{c\ 1 c\ l C2 l C2 ] +H.c). (2) 



and H e i- p h is the adiabatic Hamiltonian introduced in Ref.10, which generalizes 
the SSH Hamiltonian. 

H ep = g[ 12) ^(n Xff + n 2a )(u2 - u x ) + 7 < 12 > £>L c 2a + cL c iJ(«2 - (3) 

a a 

For symmetry reasons, the deformation ui on site i at = [(— l l )a/2, 0, 0] (i = 1, 2) 
has to be along the a:-axis, therefore U2 — u\ is the variation of the dimer length 
a, while U\ + u 2 — 0. By assuming a single phonon frequency Q, quantization is 

performed as usual by writing ui = L (b\ + b^j with L = y/h/2Am. It is convenient 

to introduce the odd-parity phonon operators d, & according to b± — b)> = \/2d^ . 
The even-parity operators do not couple to the electrons, and contribute to the 
energy a constant term. Then, after rescaling the coupling constants by 7ifi/v2£, 
H ep reads: 



Hep — 



,(12) 



X>1" + + 7 (12) J2(°l^2a + <L<ho) 



(4) 



while free oscillator term reads H p h = Ml (cftd + 1/2) To obtain an effective electron- 
only Hamiltonian, describing fermions dressed by phonons, we perform on H e i + 
H ep + Hph a unitary " displacement " transformation exp (6R) , where (from now on 

we drop the site indexes on g^ 12 ^ and 712 for short): 



R = 



g{ni +n 2 ) +l(c\ 



Co + C 9 Ci 



(d f -d) = -i?t 



(5) 
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with 8 an up to now undetermined parameter. The "displaced" Bose operators 



D 



(t) 



d(t) 



e SR are: 



"-2 J 



) +7 (c\ c 2 + 4 C 1 



(6) 



For the "displaced" Fermi operators fj = e SRs c\e SRs we have that the equations 
of motion (EOM): 



9fl/d8=(gfl +1 fl\ (dt-tf) 
dti/d8=(gfl+7fl) (dt-d) 



(7) 



can be decoupled due to the finiteness of the system, and an analytical solution is 
obtained. Taking into account the boundary conditions for 8 — 0, defining for short 
B = 8(<P — d) = — B* and indicating by Ch(x) and Sh (x) the hyperbolic cosine 
and sine, the solutions can be written: 



f! = exp(gB) \c\Ch ( 7 B) + c+Sh ( 7 B) 



/ 2 f = exp( ff B) [4Sh ( 7 5) + c^Ch ( 7 B)] (8) 

One also has ^ = ^//^ (i = 1,2). Applying Eqs.|| and || to _ff e ; + H ep + H p h 

yields the "displaced" Hamiltonian, still containing both Fermi and Bose operators. 
At this point we make the approximation of factorizing the true wave function 
\^f) into the product of Fermi and Bose functions, i.e. |^) = The 
Bose operators are now eliminated from the "displaced" Hamiltonian by averaging 
e 5R He~ 5R over a squeezed phonon wavefunction \^b) = exp(— S)\O p h) with S = 
a{d) d) — dd) and d\O p h) = 0. In the resulting one-particle polaronic Hamiltonian 
H* = (0 p h\e s e SR H e~ 5R e~ s \ p h) the renormalized interactions read: 

e* = e - hnS{2 - 8)(g 2 + ~/ 2 ) 
t* = t- 2MW(2 - <% 7 

U* = £7(3 + t 4 )/4 + (W - J)(l - t 4 )/16 - 2fift<S(2 - 8)g 2 

V* = V(7 + r 4 )/8 + (4C7 + 3J)(1 - r 4 )/32 - 2«W(2 - 6)(g 2 - 7 2 ) 

J 7 * = J(5 + 3r 4 ) /8 - (U - V) (1 - t 4 )/2 + 4MM(2 - 6) 7 2 iy j 

J*, = J(7 + r 4 )/8 - (£/ - V)(l - r 4 )/2 - 2%S18{2 - % 2 

P* = P(9 - r 4 )/8 - (C/ - V)(l - r 4 )/4 - 2MM(2 - £) 7 2 

AT* = tI - 2MW(2 - <5).g 7 



where we have defined r = (0 p ? l |e' s Ch(G'B)|e _ ' s |0pft,) = exp[— 2 7 2 <5 2 exp(— 4a)], and 
we distinguish the longitudinal (J*) and transverse (J^ y ) magnetic couplings. Our 
results differ qualitatively under several aspects from those following quasi-classical 
or perturbative treatments ~ Pj ' || ~ f§] ■ First, we find that not only the hopping 
amplitude t, but all the interactions, are non-linearly renormalized. Second, if only 
7 is assumed as non- vanishing Q then t* — t i.e. in the dimer the SSH interaction 
alone does not modify at all the hopping amplitude. Third, due to 0. < 8 < 1., the 
cooperative action of both g and 7 causes t* < t(< 0.). i.e. it promotes the charge 
carriers itineracy. 
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Figure 1: The bare (open dots) and renormalized (asterisks) electronic interactions 
for the case of two particles. 
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a (A) 

Figure 2: The electron- phonon interactions g (full line) and 7 (dotted line) vs. a. 

3 Numerical results. 

To study the ground state of the system as the dimcr length a is varied, we shall 
adopt the electron interaction parameters e, t, X, U, V, J z = J xy = P as evaluated 
as in Ref.ll, by assuming a non-degenerate orbital described by Wannier functions 
built from atomic orbitals of Gaussian shape. Also the variation with a of g and 
7 has been taken into account, following Ref.10. The effective Hamiltonian can 
be exactly diagonalized pd[ y ielding the eigenvalues as functions of the variational 
parameters 5, a and the Wannier functions shape-defining parameter pd[ T. To 
evaluate the energy of the system, one has to add to ($|_ff*|$), the contribution of 
the squeezed phonons: E p h = ?if2[Sh(2o!) 2 -|-l/2]. From the variational optimization 
of the total (i.e. electronic plus phononic) energy we obtain the true ground state. 
For lack of space we discuss only the case of two particles in the orbital. 

The bare (open dots) and effective (asterisks) interactions are shown in Fig.l 
for a phonon frequency Ml = O.leV. One notices that only at low a there is a 
strong renormalisation in the effective interactions. It is due to the fact that, as 
shown in Fig. 2, the electron-phonon interactions decrease with a, so that, at low 
a their renormalizing effect is strong, being weaker for larger a. The ground state, 
in the notation of Ref.ll, is the singlet state \Sb) for a < 2.16A, and the 5 = 
component of the triplet \T, 0) for larger a. At the transition value of a, also r and 
r vary, indicating that the system interactively readjusts both the phononic and 
the electronic features. 

Knowing the eigenstates |$) of H* we can evaluate the correlation function (CF) 

((Y)) for an operator Y according to ({Y}) = (<I>\(0 ph \e s e SR Ye- SR e- s \0 ph )\<S>) . By 
defining 9 as in Ref.ll, one has in the \Sb) state: 

((^1^2)) = cos 2 9 

«K T n u + n 2] n 2l ))) = {{{^c ] n c 2l c 2] + H.c.))) = [1 - r 4 cos(20)]/2 (10) 
(((S+S 2 - + S^S+)))=mSfSi)) = -[l + r 4 cos(20)]/2 

Fig. 3 shows the a dependence of the charge distribution CF ((77,1^2)), of the trans- 
verse magnetic CF (((Si S 2 + H.c.))) , from which the longitudinal magnetic CF can 
be deduced by scaling, and of the on-site bipolaron (OSB) CF (((n 1 ^n 1 ^ +n 2 ^n 2 ^))) 



5 



'1 1.5 2 2.5 3 3.5 

a (A) 

Figure 3: Correlation functions. See the text for explanation. 



coinciding with the OSB hopping CF {{{c\^c\^c 2 yC 2 ^ +H.c.))). For a small, where g 
and 7 have a strong effect, the formation of an OSB is possible, implying an uneven 
distribution of the charges between the sites, and, at the same time, the magnetic 
CF's are sizeable, even though not fully developed. The phonons therefore induce 
the cohexistence of charge and magnetic correlations for dimer lengths such that 
the electron-phonon interactions are not negligible, and promote the itineracy of 
the OSB. 
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